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Abstract 

We study the effects of twisted boundary conditions on the quark fields in the epsilon regime of 
chiral perturbation theory. We consider the SU{2)l x SU(2)r chiral theory with non-degenerate 
quarks and the SU(3)l X SU(3)r chiral theory with massless up and down quarks and massive 
strange quarks. The partition function and condensate are derived for each theory. Because flavor- 
neutral Goldstone bosons are unaffected by twisted boundary conditions chiral symmetry is still 
restored in finite volumes. The dependence of the condensate on the twisting parameters can 
be used to extract the pion decay constant from simulations in the epsilon regime. The relative 
contribution to the partition function from sectors of different topological charge is numerically 
insensitive to twisted boundary conditions. 
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I. INTRODUCTION 



Understanding the strongly interacting regime of Quantum Chromodynamics (QCD) re- 
mains a challenging problem in physics. The dynamics is non-perturbative and results in 
the confinement of the colored quark and gluon degrees of freedom. Numerically simulating 
the gauge theory on a spacetime lattice provides a first principles method for determining 
the properties of the color neutral hadrons that are observed in nature. Recently, there has 
been much effort in applying effective field theory methods to quantitatively compute the 
effect of lattice artifacts such as unphysical quark masses, quenching, finite lattice spacing 
and finite volume on lattice simulations of hadronic observables. Effective field theory cal- 
culations can be used to extrapolate the results of simulations with unphysical parameters 
to real world observables |l( . One can also use the dependence on lattice artifacts to extract 
the parameters appearing in low energy effective Lagrangians of QCD. 

Lattice QCD simulations usually employ periodic boundary conditions for the quark and 
gluon fields. The results of a simulation should be independent of the choice of boundary 
conditions as one takes the volume to infinity, so there is considerable freedom in choosing 
the boundary conditions for the fields. However, the action must be single valued so that 
observables are well-defined. Fields need only be periodic up to a transformation which is 
symmetry of the action. If a generic field <fi depends on a lattice coordinate x% which is 
defined on the interval < Xi < L, then we must have 

<f)(xi + L) = U<p(xi) , 

where U is a symmetry of the action and U'U — 1. When U ^ 1, these are called twisted 
boundary conditions. This is equivalent to gauging a global symmetry of the action and then 
turning on a uniform background gauge field. In infinite space, a constant background gauge 
field can be removed by a gauge transformation, so there are no physical consequences in the 
infinite volume limit. However, a uniform background gauge field can modify the physics 
of a compact space. Different applications of twisted boundary conditions have appeared 
throughout the years 0, 0, HL M, IE EL El Hi- 1 Recently, there has been a renewed 



interest in twisted boundary conditions in lattice QCD simulations 
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recent work has focused on using twisted boundary conditions on quarks as a means of 
avoiding the restriction of lattice momenta to integer multiples of 2n/L. 

In this paper, we consider a different aspect of twisted boundary conditions. Since dy- 
namical effects due to twisting show up as corrections at finite volume, it is natural to 
examine the effects of twisted boundary conditions in the ^-regime of chiral perturbation 
theory (xPT), where volume effects are non-negligible. ;^PT is the low-energy effective 
theory of QCD that encodes the dynamics of pseudo-Goldstone bosons that emerge from 
spontaneous breaking of chiral symmetry. Spontaneous symmetry breaking does not occur 
in finite volume. Therefore the perturbative expansion of x?T must breakdown in the limit 
m q — > with L held constant. Indeed when m^L <C 1 the zero momentum mode of the 
Goldstone boson becomes strongly coupled. In this limit, the path integral over the zero 
modes must be treated exactly, while nonzero modes can be treated perturbatively. Com- 
putation of the chiral symmetry breaking order parameter in this regime shows that the 



1 The use of twisted boundary conditions in the time direction can be utilized to produce an imaginary 
chemical potential 
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chiral condensate vanishes exponentially as m%.f%(3V — > 0, where f3V is the four- dimensional 
volume of the lattice. For early references on these finite-size effects in theories with spon- 
taneously broken symmetries, see |25l . l26l. 271. 28l. 129. I30L l31l . For later investigations and 




numerical simulations see, for example, |32|, |33j, |3J, |35|, |3JJ |37|, |3J 

In the p-regime of chiral perturbation theory one takes 1/L ~ m v ~ Q so propagators 
scale as Q~ 2 , vertices scale as Q 2 or higher (since vertices have at least two derivatives or 
one insertion of the quark mass) and loop integrations (or sums in finite volume) scale as 
Q 4 . The contribution to the partition function from an arbitrary vacuum graph scales as 
Q n where n > 41 — 21 + 2V = 2/ + 2, / is the number of loops, I is the number of internal 
lines and V is the number of vertices. Thus, in this regime the partition function can be 
computed in perturbation theory. In the e-regime, 1/L ~ e but m n ~ e 2 . Thus, the zero 
momentum mode is singled out because its propagator scales as e~ 2 while the nonzero mode 
propagators scale as e -1 . The leading order interactions of zero modes with themselves are 
proportional to m 2 ~ e 4 so an arbitrary vacuum graph with only zero momentum modes 
and leading order vertices scales as e 4l - 41+4V = e 4 . Therefore, an infinite number of graphs 
with zero momentum modes must be resummed at leading order. This resummation is 
straightforward because the theory of zero modes alone is simply a matrix integral. 

Now we see how this power counting is modified in the presence of twisted boundary 
conditions. The twisted boundary conditions induce a constant shift in the momentum of 
the propagators of fields that are charged under the background gauge field: 

(1) 



k 2 + ml (k^ + B^) 2 + ml ' 

where is the background gauge field. For generic twists, B^ ~ 1/L ~ e and the propaga- 
tor in Eq. (0) scales as e~ x . The zero modes of fields charged under the background gauge 
field are no longer strongly coupled. However, if there are Goldstone bosons that are not 
charged under the background gauge field their zero modes will continue to be strongly cou- 
pled. In QCD with N quarks, B^ must couple either to baryon number j20|, in which case 
all meson fields are neutral, or an element of the SU(N) flavor group. Because quark mass 
differences and electric charge break the flavor symmetry down to the diagonal subgroup, 
we will only consider twists that are flavor diagonal. Then the flavor-neutral Goldstone 
bosons will be neutral under the background field and there will still be a restoration of 
chiral symmetry as m q — > with L fixed. 2 The qualitative behavior of the condensate as 
a function of m q is the same as without twisting. However, not all the zero modes of the 
coset SU(N)l x SU(N)r/SU(N)v are strongly coupled, only those valued in the diagonal 
subgroup. The behavior of the condensate is quantitatively different in this case. For suf- 
ficiently small twists the background gauge field can be comparable in size to m 2 and the 
appropriate power counting is B^ ~ m n ~ e 2 . In this case the zero modes of all Goldstone 
bosons are strongly coupled. We will refer to the case of B^ ~ m n ~ e 2 as the twisted 
£-regime and B^ ~ e, m n ~ e 2 as the remnant e-regime. 

In the untwisted theory at lowest order in the e expansion, the partition function depends 
only on the mass term in the chiral Lagrangian. Computation of the volume dependence of 
the chiral condensate can be used to extract m q (qq) OQ , where (qq)oo is the infinite volume 



2 Even for generic twists it is easy to show that there will be N — 1 Goldstone bosons neutral under the 
background gauge field. 
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condensate. In the presence of twisted boundary conditions, however, the kinetic term 
in the chiral Lagrangian gives a leading order contribution to the partition function. By 
computing the finite volume condensate as a function of both volume and twist, it is in 
principle possible to extract both rriq^qq)^ and /„-, where f n is the pion decay constant. 
This latter quantity is usually extracted from the late time behavior of the axial correlator. 
Thus, twisted boundary conditions give an alternative method by which one may extract f n 
from simulations in the e-regime. 

We also evaluate the projection of the partition function onto sectors with fixed topolog- 
ical charge by including a #-term in the chiral Lagrangian and Fourier transforming with 
respect to the parameter. This calculation can be done analytically for the untwisted 
theory j42| . We perform the calculation for the isospin symmetric SU{2)l x 577(2) r chiral 
theory. Twisted boundary conditions do not dramatically change the results of Ref . |42j . In 
particular, the observation that the QCD partition function is dominated by the sector with 
zero topological charge when m g (qq) oo pV < 1 still holds in the presence of twisted boundary 
conditions. 

The paper is organized as follows: in SectionlHl we describe in detail the twisted boundary 
conditions imposed on the quark fields in QCD with N quark flavors and the twisted finite- 
volume meson propagator. We write down a general formula for the partition function 
in the presence of twisted boundary conditions. In Sections IIHI and IIVI we evaluate this 
partition function as well as the finite volume condensate for twisted SU(2)l x SU(2)r and 
SU(3)l x 577(3) # chiral theories, respectively. We consider both the twisted e-regime and 
the remnant e-regime. In Section |3 we compute the projection of the partition function 
of the 577(2) £ x SU(2) R chiral theory onto sectors of fixed topological charge. Section IVD 
contains a brief summary and discussion of directions for future work. 



II. GENERALITIES 

Let us consider lattice QCD gauge theory with N quark flavors. In a finite volume, 
the boundary conditions on the quark fields must be chosen so that action is single-valued. 
Thus the generator of any flavor symmetry of the action can be used to specify the boundary 
conditions of the quark fields. In the i th spatial direction, we impose flavor diagonal twisted 
boundary conditions of the form 

q(x l + L) = e^q(x l ), (2) 

where the Abelian phases M can be expanded in the Cartan subalgebra of U(N). Instead 
of using the Cartan generator basis, we will sometimes choose to write the decomposition of 
twists in the flavor basis: (Q^ab = 8^9°, where 9% = for all flavors. Furthermore due to 
periodicity, we need only consider twist angles Of in the range — n < Of < n. 
It is convenient to introduce quark fields q(x) defined by 

q(x) = W(x)q(x), (3) 

with 

U(x) = e ie ^ L . (4) 

The q(x) fields satisfy periodic boundary conditions. The quark part of the QCD Lagrangian 
can then be rewritten in the form 

£ = j(fl+M)q, (5) 
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where the effect of twisting appears as a background gauge field 



D^ = d^ + iBp, (6) 

where = Q^/L. 

To describe QCD at low energies, one turns to an effective field theory, chiral perturbation 
theory, that is written in terms of the pseudo-Goldstone bosons that appear from chiral 
symmetry breaking. As shown in @, in the presence of twisting, this low energy effective 



theory of QCD can be rewritten in terms of a modified coset field given by 

£(x) = tf(x)Z(x)U(x), (7) 

that satisfies periodic boundary conditions. The coset field however, has non- vanishing 
charge in the background gauge field, B^. In terms of the modified coset field, the chiral 
Lagrangian appears as 

£ = f -f [tT(D^D^) - tr(Ex + x tt)] , (8) 
where x is proportional to the quark mass and the covariant derivative is 

D^t = d^t + i\B^t}. (9) 

From the Lagrangian in Eq. (|S|). we can deduce the propagator of the meson fields. For 
our general discussion, we begin by omitting the flavor structure of the propagator from our 
consideration. In a finite volume with dimensions L x L x L x (3, the meson propagator for 
the mode = (n, 714) is proportional to 

(10) 



(27m 4 //3) 2 + (27rn + 0) 2 / L 2 + ml ' 



where n G Z 3 and n 4 £Z due to the periodic boundary conditions and the term arises from 
the background gauge field. For power counting we will always assume (3 ~ 1/L ~ e. When 
one works in the e-regime, where m^L 1, the zero mode (n^ = 0) of the finite- volume 
propagator ordinarily is strongly coupled. Twisting however has introduced an effective 
mass squared B 2 + m 2 for the zero mode. Because the flavor-neutral mesons are unaffected 
by the twisting, there are always strongly coupled zero modes in the group manifold. In the 
twisted e-regime all Goldstone boson zero modes are strongly coupled. Since the zero modes 
must be treated non-perturbatively, we separate the field into two pieces 

t(x) = t e 2i ^ x),f \ (11) 

where So is the zero mode and the <fi contains the non-zero modes of the pion fields. 

Returning to the Lagrangian in Eq. (jHJ), we see that to lowest order the partition function 
in the twisted e-regime is 

Z = Jvt exp j^Str^So^So + Ke( X S )]} , (12) 



where V — L . We have dropped a B^ dependent constant which only modifies the nor- 
malization of the partition function and is physically irrelevant. Note that this partition 
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function is very similar in form to that which arises when considering QCD with an isospin 
chemical potential 0, 44[ 45|. However, an isospin chemical potential induces a that is 



time-like rather than space-like and the contribution to the action comes with the opposite 
sign. For our treatment of the zero modes in the SU{2) L x SU{2) R and SU(3)l x SU(3)r 
chiral theories, the expression for the partition function in Eq. (|12|) will be our starting 
point. We will also consider the partition function for the remnant e- regime, where only the 
flavor-neutral mesons have strongly coupled zero modes. In this case, the matrix E belongs 
to the diagonal sub- manifold and only the quark mass term survives in Eq. (JT2J). 



III. TWISTED e-REGIME OF SU(2) 

In this section we consider the effects of twisted boundary conditions on the quark fields 
in two flavor xPT in the twisted e-regime. The simplicity of the group manifold will allow us 
to reduce the zero mode contribution to the partition function to a one-dimensional integral 
which is easily evaluated numerically. The dependence of the finite- volume quark condensate 
on twisting angle is then investigated. We also derive the partition function in the remnant 
e regime when only the ir° zero mode is strongly coupled. 

For two quark flavors, the general form of the background gauge field is 

b, = diag(5;, b% = \(b; + b$ + \(b;- b'k (is) 

The matrix integral appearing in Eq. (jl2j) can largely be carried out due to the simple 
structure of the 577(2) group manifold. First notice the term proportional to the identity in 
Eq. (fH?j) decouples from the zero mode integration. To perform the remaining integration 
in the twisted e-regime, we use a familiar parametrization of SU (2) 

S = cos a + in ■ crsina, (14) 

where h is a three-dimensional unit vector. With this choice for the arbitrary element of 
577(2), the group measure, normalized to unity, is 

£>Eo = — g / dQn I da sin 2 a, (15) 

and the partition function is straightforwardly evaluated. 

The final result can be cast into the form of a one- dimensional integral over well known 
functions. To this end, we define the dimensionless variables s and t 

s = \f^Vtx{ X ) = \f 2 JVml = -^pt(qq)ooPV, (16) 

t = \{B U - B d ffjV = {6U ~ L f )2 f 2 JV, (17) 

where m v is the pion mass in the infinite volume limit. The variable s is the ordinary 
parameter relevant to the e-regime 27j ]. while the variable t contains the new effects due to 



twisting the quark fields at the boundary of the lattice. Notice t is zero only if in each spatial 
direction the quark phases at the boundary are pairwise identical. Furthermore in evaluating 
the matrix integral, we do not need to assume that the mass matrix is proportional to the 
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identity. All results at this order in the chiral expansion depend only on tr(%). Terms in 
lead to a departure from this result, but these are beyond the order we work. 
The angular integral in the group measure can be done analytically. In terms of the 
variables t and s, the partition function to leading order in the twisted e-regime is 

1 f 7 ^ sin ot 

Z(s,t) = —= / da — -=- Erfi(\/t sin a) exp(2 s cos a — t sin 2 a), (18) 
V 71 " Jo Vt 

where Erfi(x) = Eri(ix)/i, and Erf(x) is the standard error function. Untwisting the quark 
fields, i.e. for t = 0, we arrive at 

Z(s, 0) = — / da exp(2 s cos a) sin 2 a, (19) 
TT Jo 



which is the form of the ordinary SU(2) group integral in the e-regime |27j. In terms of a 
modified Bessel function, we have Z(s, 0) = I\{2s)/s. 

In the twisted e-regime, ~ m n ~ e 2 and t/s ~ 1, but in the remnant e-regime 



s mi 



So the large t limit corresponds to the remnant e-regime. Taking the asymptotic limit of 
the error function, 

Erfiv/x e x , . 

hm = — , (21) 



X 



we find the t-dependence in Eq. ()18|) reduces to an overall multiplicative factor which can 
be pulled outside the integral. In general, one is interested in computing expectation val- 
ues of operators which can be calculated by adding sources to the Lagrangian and then 
differentiating the logarithm of the partition function with respect to these sources. The 
condensate is obtained by differentiating log Z(s,t) with respect to s (see Eq. (|2^j) below). 
So the multiplicative dependence on t is a trivial modification of the normalization of the 
path integral which can be ignored. Thus, for large t, the partition function is 

Z(s) oc / da exp(2 s cosa) oc Io(2s) . (22) 
Jo 

The 7T° zero mode has a mass ~ m n ~ e 2 but the charged pions have an effective mass 
~ \B^\ ro e. Physically, one expects to be able to integrate out the charged pion zero modes 
in this limit. The twisted e- regime partition function should reduce to that of the remnant 
e-regime, in which only the 7r° zero mode is dynamical. The remnant e-regime partition 
function is obtained using the parametrization 

S = exp(zacr 3 ), (23) 

and the group integral measure is trivial since it is an integral over a U(l) subgroup of the 
coset SU(2)l x SU{2)r/ 'SUv(2). This immediately leads to a partition function identical 
to that in Eq. 
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FIG. 1: The quark condensate of the SU(2)i x SU(2)r theory with twisted boundary conditions 
is plotted as a function of s, for several values of the twisting parameter t. Results in the remnant 
e-regime correspond to the curve denoted t = oo. 



The quark condensates for the up and down quarks are identical, (uu) = (dd), because 
the partition function depends only on the trace of the mass matrix tr(%). We denote either 
of the condensates by (qq). They are calculated by differentiating the partition function 

SHI*™ < M » 

In Figure ^ we show the dependence of the quark condensate on the twisting parameter t. 
For t — 0, the result is just that of the ordinary e-regime of SU{2). The condensate decreases 
with s, and vanishes as s —>■ due to the restoration of chiral symmetry at finite volume. 
For t ^ there are twisted boundary conditions on the quark fields and the qualitative 
behavior of the condensate as a function of s is unchanged. The value of the finite-volume 
condensate, however, is altered and depends on the twisting variable t. As t becomes larger, 
the effective mass of the charged pion zero modes increases due to the twisting and these 
modes can no longer conspire to restore chiral symmetry. Therefore the condensate increases 
in magnitude with t for fixed s. When t/s is very large, the theory is in the remnant e-regime 
and only the ir° zero mode is relevant. The result for the remnant e-regime is denoted by 
the curve t = oo. We see that the behavior of the condensate for non-zero t approaches 
that of the remnant e-regime in the limit of large t. In Figure 121 we plot the dependence 
of the quark condensate as a function of t for fixed s. Mapping out the dependence of the 
finite-volume condensate as a function of U — 6 d for fixed four-volume and s gives one an 
alternate method to determine /„.. As s — > the condensate must vanish because of the 
restoration of chiral symmetry and therefore the value of the condensate is insensitive to 
the twisted boundary conditions. As s — > oo one is approaching the infinite volume limit 
and all observables should be insensitive to boundary conditions. The condensate is most 
sensitive to twisted boundary conditions near s — 1. 
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FIG. 2: The quark condensate of the SU(2)l x SU(2)r theory with twisted boundary conditions 
is plotted as a function of the twisting parameter t, for several values of s. 



IV. TWISTED e-REGIME OF 517(3) 

In this section, we investigate the consequences of twisted boundary conditions on QCD 
with three flavors of quark fields. The group integration over the ST/ (3) manifold is more 
complicated than 577(2) . Our formulae for the matrix integral in the most general case is an 
eight dimensional integral which we do not know how to perform analytically. To simplify 
numerical analysis of the partition function, we work in the limit of massless up and down 
quarks and choose isospin symmetric boundary conditions. The partition function in this 
case can be reduced to a one-dimensional integral and the dependence of the strangeness 
condensate on twisting angles and strange quark mass can be easily calculated. 

For the case of the 577(3) flavor group, the flavor- diagonal background gauge field is 

B, = diag(fl£, B% B$ = Bl + 5jA 3 + SjA 8 , (25) 

where 5° = §(B£ + B* + S'), B* = §(£« - B% and Sj = ^(5« + - 2S'). In 
Eq. (|12|). any terms involving 5° decouple from the group integral. For the arbitrary group 
element of SU(3), we use a convenient form jlril . 47] that is a generalization of the Euler 
angle decomposition, 

With this decomposition, the invariant volume element of the S77(3) group is ji^, |4?i IH[ 

PE oc / / dad"fdadc / dfidbdip sin 2/? sin 26 sin 2^ sin 2 ?/' . (27) 

Jo Jo Jo 
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(Note that in this case we have not chosen to normalize the group integral to 1.) 

For the moment, let us consider three independent flavor twists and three non-degenerate 
quark masses. In this general case, we can write the partition function in the twisted e- 
regime as a function of the quark mass dependent variables s a = \x a f 2 p3V and the various 
combinations of flavor twists, t ab = \B a ■ B b f 2 (3V. We find the partition function has the 
form 

Z(s a ,t ab ) = N I Vt exp[f(s u ,s d ,s s ) + g(t 33 ,t 38 ,t 88 )], (28) 



where the quark mass dependent term in Eq. (fT2j) contributes to the function f(s u , s d , s s ), 
which is given by 



r I U d 

f(s ,s ,s 



+ 



s COS ip COS — p= 



s cos a + c + a + 7 — 



s cos a — c — a + 7 



+ s u cos ip cos (a + c + a + 7 + 



+ s cos ip cos (a — c — a + 7 + 



(29) 
cos b cos (3 
sin b sin/?, 



and the term in Eq. (fT^) which depends on the background gauge field contributes to the 
function g(t 33 ,t 38 ,t 88 ), which is given by 



g(t 33 ,t i8 ,t 88 ) = -t 66 [(3 + cos 2i/j) cos 2(3 cos 26 — 4 cos(2a + 27) cos ip sin 2(3 sin 2b] 

3 

- V3 t 38 (cos 2(3 + cos 2b) sin 2 ip + - t 88 cos 2ip. 



(30) 



The full expression for the partition function in Eq. (|28j) is an eight dimensional integral 
which appears difficult to evaluate analytically in the most general case. We will make two 
assumptions which greatly simplify the evaluation of the partition function. It is convenient 
to choose isospin symmetric boundary conditions with B d = B™. In this case, we have 



t' 

git') = 2 cos2 ^' 



(31) 



where t' = \{B U — B s ) 2 f%/3V . Secondly, we make the physically relevant approximation 
of neglecting up and down quarks masses while letting the strange quark mass be nonzero. 
Then the action only depends on the angles (p and ip and all but one integral can be done 
analytically: 



Z(s s ,t' 



N' 



V3n 



w/2 



dip sin 2ip sm z ip exp I s s cos ip cos —= H — cos 2ip 



N" / dxx{l-x 2 )I {s s x)e t,{x2 - l/2 \ 



(32) 



where Iq{x) is a modified Bessel function. 

The strangeness condensate, (ss), is given by 



(ss) 
(ss), 



_d_ 

dp 



logZ(s s ,t'), 



(33) 
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FIG. 3: The strangeness condensate in the SU(3)l x SU(3)r theory with twisted boundary con- 
ditions as a function of s s , for several values of the twisting parameter if . The remnant e-regime 
corresponds to the curve denoted if = oo. Note that B u = Bd and m u = = 0. 

where the corresponding condensate in the infinite volume limit is denoted by (ss)oo- I n 
Figure El we show the dependence of the strangeness condensate on the twisting parameter 
if . The result for t' = 0, corresponding to the case of periodic boundary conditions, is the 
finite-volume behavior of the condensate in the ordinary e-regime. Like the (qq) condensate 
in the SU (2) L x SU (2) R theory, the (ss) condensate decreases monotonically with m s . When 
t' is non-zero, the theory is in the twisted e-regime and the qualitative behavior as a function 
of s s is unchanged. However, the flavor non-diagonal zero modes become less important than 
the flavor diagonal modes. The condensate is larger in the presence of twisted boundary 
conditions. The greatest enhancement is achieved when t'/s s ^> 1, which is the remnant 
e-regime. 

Although we have not been able to analytically derive the large if limit of the strangeness 
condensate directly from Eqs. (}3*2*|) and ()33|) . we can treat the remnant e-regime exactly 
when m u = rrid = 0. In the remnant e-regime, only the 7T° and rj zero modes remain strongly 
coupled. The integral over the SU(3)l x SU(3)r/ SU(3)v coset reduces to integrals over 
the U(l) diagonal subgroups. Therefore, the measure is trivial and S can be parametrized 
by the diagonal matrix S = di&g(e ldl , e 1 ® 2 , e~ l( - ei+e ^) . Thus, the partition function in the 
remnant e-regime is 



To compare with our numerical results in the twisted e-regime, we set s u = s d = 0, and find 
Z(s s ) oc Iq(s s ). The resulting strangeness condensate is labeled by if = oo in Figure El 



Z(s a ) — N f 




s u cos 9 1 + s d cos 9 2 + s s cos(#i + 2 )] 



(34) 
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V. TWISTED QUARKS IN FIXED TOPOLOGY 



In this section, we discuss the effect of twisted boundary conditions on the partition 
function in sectors of fixed topological charge. If Z v corresponds to the partition function in a 
sector with topological charge v then Z v jZ can be interpreted as the probability of finding the 
gauge configuration with topological charge v. This can be computed by introducing a 9 term 
in the chiral Lagrangian and computing the Fourier transform of the partition function with 
respect to 6. Ref. 42| showed that the QCD partition function is dominated by vanishing 



or small v in the e-regime. We study whether this conclusion is modified in the presence of 
twisted boundary conditions. For simplicity we consider only the SU{2)l x SU{2)r theory 
with m„ = m<f. 

In the twisted e-regime, the partition function for vacuum angle 9 is given by 

Z{6) = j Vt exp j^tr^Eo^Ej + $te(e^ 2 X t )} } , (35) 

where we have specialized to the case of two flavors. In a sector with fixed topological charge 
u, the corresponding partition function is 

*^ e -**Z(0), (36) 

and can be straightforwardly evaluated. We find 

Z u (s,t) = ^ da ^3 Erfi( v / t sin a) e~ tsin2 a I 2u {2s cos a). (37) 
V 71 " Jo vt 

Note that Z u (s,t) = Z^ u (s,t). It is easy to show that in the limit t — > 0, we recover the 
result 

Z v (s, 0) = / I/ (s) 2 -/ 1/+1 (s)/ I/ _i(s). 

In Figure|3]the logarithm of the ratio Z u (s,t)/Z(s,t) is plotted for v = through v = 4. In 
this plot we show the result for two different values of s with and without twisted boundary 
conditions. We see that the dominance of the sector with zero topological charge in the 
e-regime is not modified by twisting. The numerical value of Z u (s,t)/Z(s,t) changes little 
between t — and t = 10. 



VI. SUMMARY 

In this paper, we have investigated the effects of twisted boundary conditions on the quark 
condensate in small volumes. We worked in the twisted e-regime of chiral perturbation 
theory where 1/L ~ e, ~ m n ~ e 2 and all Goldstone boson zero modes are strongly 
coupled. The partition function and quark condensate were computed to leading order by 
performing the path integral over these zero modes. We also investigated the remnant e- 
regime, where 1/L ~ ~ e, m n ~ e 2 . In this case only the flavor-neutral mesons are 
strongly coupled and the contribution from these modes dominates the partition function. 
It is the flavor-neutral zero modes that lead to a restoration of chiral symmetry at small 
volumes in the remnant e-regime. 
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FIG. 4: The partition function for fixed topological charge. We plot logi [Z U (s,t)/Z(s,t)] as a 
function of v for different s and t. The solid lines are provided merely to guide the eye. 



The value of the quark condensate in small volumes affects the masses of particles. For 
example, the pion mass is given by 



in 



2 



tr(x) (qq) 
2 (qq)o 



(38) 



and the masses of heavy particles such as the nucleon |49j , pick up an additive renormaliza- 
tion, 

m N = m -a tr(x) tPt— • (39) 

W)oc 

When one imposes twisted boundary conditions on the quark fields, particle masses are 
modified since (qq) is larger relative to its value in the absence of twisting. Of course, the 
pion mass vanishes at very small volumes in any case. 

Lattice calculation of the dependence of quark condensates and particle masses on twisted 
boundary conditions can be used to test the predictions of twisted chiral perturbation theory. 
Moreover because the kinetic energy term contributes to the twisted partition function, it 
is in principle possible to determine both the infinite volume condensate and the pion decay 
constant from numerically determining the finite volume condensate as a function of the 
twist angles and volume. Lastly, we showed that sectors with zero topological charge remain 
the dominant contributors to the partition function in the twisted and remnant e-regimes. 

For the results in this paper to be useful in actual simulations, it may be necessary to 
extend our calculations to include effects due to quenching or partial quenching as well as 
perturbative loop corrections from nonzero modes. Note that the simplicity of the zero-mode 
manifold in the remnant e-regime may make it easier to compute the partition function for 
quenched and partially quenched theories with twisted quarks. It would also be interesting 
to examine other correlation functions, especially two point functions of currents. We hope 
this paper will stimulate further numerical and analytical work in this novel regime of QCD. 
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